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UNSTABLE GRB PHOTOSPHERES AND E ± ANNIHILATION LINES 
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ABSTRACT 

We propose an emission mechanism of prompt gamma-ray bursts (GRBs) that can reproduce the observed 
non-thermal spectra with high radiative efficiencies, > 50%. Internal dissipation below a photosphere can 
create a radiation-dominated thermal fireball. If e ± pairs outnumber protons, radiative acceleration of e ± pairs 
drives the two-stream instabilities between pairs and protons, leading to the "proton sedimentation" in the 
accelerating pair frame. Pairs are continuously shock heated by proton clumps, scattering the thermal photons 
into the broken power-law shape, with the non-thermal energy that is comparable to the proton kinetic energy, 
consistent with observations. Pair photospheres become unstable around the radius of the progenitor star where 
strong thermalization occurs, if parameters satisfy the observed spectral (Yonetoku) relation. Pair annihilation 
lines are predicted above continua, which could be verified by GLAST. 

Subject headings: gamma rays: bursts — gamma rays: theory — radiation mechanism: non-thermal 



1. INTRODUCTION AND SUMMARY 

The emission mechanism of the prompt Gamma-ray bursts 
(GRBs) (the most luminous objects in the universe) is 
still enigmatic despit e the recent pr ogresses in the Swift 
era (iMeszarosI l2006t IZhang 120071). Main issues are 
(A) t h e efficiency probl em dlokaetaT] 12006: Zhan g et ail 
12001: iTomaetaLl [20061) and (Bi the cooling problem 
(iGhisellini. Celotti & Lazzatill2000t IMeszaros & Reesll2000l) . 
The first problem (A) is that the GRB radiative efficiency, de- 
fined by the GRB energy divided by the total energy includ- 
ing the afterglow energy, is too high (> 50%) to be produced 
by the internal shocks. Although the high efficiency may be 
achieved b y a large dispersion in the Lorentz factor of the 
outflows dBeloborodovl 120001: iKobayashi & Sarill200lh . this 
would yield smaller correlation coefficient for the observed 
spectral correlations, such as Ama ti and Yonetoku relations 
(lAmatill2006H Yonetoku et al.l l2004). as far as the non-thermal 
(synchrotron or inverse Compton (IC)) processes determine 
the spectral peak energy. The recent Swift observations make 
the problem even worse since the early afterglow energy is 
smaller than expected, raising some of GRB efficiencies up 
to > 90%. The second problem (B) is that the cooling time is 
much shorter than the dynamical time, making the low-energy 
spectral slope steeper than the observations. 

Motivated by these problems, thermal photosphere models 
are proposed where the outflow energy is in ternally dissipated 
and thermalized inside the photosphere (|Thompsonl Il994t 
IGhisellini & Celottil 119991: IMeszaros & Reesl 120001) These 
models have an advantage to stabilizing t he peak energy, 
which is identified with the t hermal peak ( Thompson et al.l 
120071: iRees & Meszarosl 12005). Thermal peak s may be ac - 
tually associated with up to 30% of long GRBs dRvdell2005l) . 

However a simple photosphere model is not compatible 
with observations that most GRBs are highly non-thermal. 
This is one reason that excludes the original fireball model 
dPacz vhskil 1 9861: iGoodmanl 1 986b . Although Comptonization 
of the thermal photons via magnetic reconne ction and/or tur- 
bulence is invoked for non-ther mal spectra dThompson et al.l 
2007; Gia nnios & Spruitl 120071) . the heating mechanism of 
electrons is largely uncertain. 
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In this Letter we propose a possible scenario to make pho- 
tospheres non-thermal. We show that radiatively driven in- 
stabilities 3 occur in a radiation-dominated photosphere if e^ 
pairs outnumber protons. 4 Since radiation selectively pushes 
pairs rather than protons, the relative velocity between pairs 
and protons increases, driving the two-stream instability, pos- 
sibly of the Weibel type (§|3]l. This produces the small-scale 
inhomogeneity of the proton-to-pair ratio, that could grow up 
to proton clumps causing shocks with pairs. Shocked pairs 
scatter thermal photons into the observed non-thermal spectra 
with energy comparable to the proton kinetic energy, consis- 
tent with observations. The unstable pair photosphere may be 
comm on since the observ ations including the Yonetoku rela- 
tion (I Yonetoku et al.f 2004) suggest that the instabilities occur 
around the radius of the progenitor star where strong ther- 
malization occurs (§ |2). Our model employs the radiation- 
dominated fireball, solving (A) the efficiency problem, and 
naturally achieves the continuous heating of pairs, solving (B) 
the cooling problem (§0). Pair annihilation lines are predicted 
above continua, whi ch could enable GLA ST to verify the pair 
photosphere model dMurase & Iokall2007l) (§ EJ. We use the 
unit Icb = h = 1 and Q x = Q/\Q X in cgs units unless otherwise 
stated. 



2. PAIR PHOTOSPHERE AND YONETOKU RELATION 

Most opacity of the fireball photosphe r e can be pro- 
vided by e^- pairs (IRees & Meszarosl 120051 IMeszaros et al.l 
120021: iPe'er & Waxmanl 120041) . The observed GRB spec- 
trum, if extrapolated, has a signific ant fraction of energy 
above the pair producti on threshold (Lithw ick & Saril [20011: 
iBaring & Hard ing 1997). So more pairs can be produced than 
electrons associated with protons. (Pairs are m p /m e ~ 2000 
times more abundant if they have the same energy as pro- 
tons.) We assume that internal dissipation such as shocks 
produces pairs via non-thermal processes and create the pair 
photosphere. 

The plasma and photons are subsequently thermalized by 
scatterings under the photosphere. Assuming that we identify 
the thermal peak T t, s = TT with the peak energy of GRBs, we 
obtain the luminosity L ~ 4irr 2 aT 4 cT 2 and hence the comov- 



3 Our instability is different from that in Waxman & Piran 1 1994). 

4 We only consider protons for baryons for simplicity. 



2 



Ioka et al. 



ing size of the photosphere as 



e p = r/T~4xl0 & L l 5 { 2 T o t a cm, 



(1) 



where T b s .2 = T bs I 'lOOkeV and T is the bulk Lorentz factor. 
The optical depth of the pair photosphere is r 7C ~ n±a T l p ~ 
1 , which yields the comoving number density of pairs as 

»± -3x10%;^ 2 cm" 3 . (2) 

Our model requires that pairs are accelerated by radia- 
tion, so that the radiation dominates the pair rest energy, 
aT A > n±m e c 2 . We also require that protons are not accel- 
erated by radiation exerted on the associated electrons which 
is ~ (n p /n±)aT 4 since we are considering n± 3> n p and 
T ie ~ n±<jjl p ~ 1, so that (n p /n±)aT 4 < n p m p c 2 . With equa- 
tions (Q~|i and (O, these limit the photospheric radius as 

1 x lO'^Kf cm < r < 7 x lO'T^f cm. (3) 



If we insert the observed Yonetoku relation (Yoneto ku et al.1 
2004), 



T bs oc L 1,/2 (Yonetoku relation) 



(4) 



into the above relation, we have very weak parameter de- 
pendence as oc L" 1 / 8 . Interestingly the radius in Eq.( 3) is 
comp arable to the Wolf-Rayet stellar radius ~ 2-20,R© (ICoxl 
2000). Considering that the strong dissipation occurs within 
the progenitor star via interactions at the jet-star boundary 
(Zhang et al. 2004), most GRBs could satisfy our model as- 
sumptions. Pairs expand to the light speed within the expan- 
sion time ~ Ipjc with the comoving acceleration 



£ p /c 



2 x 10%; /2 rl i2 cm/s 2 . 



(5) 



For later use we define the mean proton density n p , which 
satisfies h p <C n± for the pair photosphere. In this Letter, we 
take n±/n p = 10 3 , T obs = 100 keV, L = 10 51 erg/s and T = 10 3 
for fiducial parameters. The energy density ratio is then aT 4 : 
n p m p c 2 : n±m e c 2 ~ 1 : 0.4 : 0.2. 

3. WEIBEL INSTABILITY BETWEEN E ± AND PROTONS 

Radiation selectively pushes pairs rather than protons be- 
cause the cross section for scattering is proportional to the 
inverse square of mass aj oc m~ 2 . Then the relative veloc- 
ity between pairs and protons must rise in the initial stage. 
Here Coulomb collisions between pairs and protons are usu- 
ally negligible. Collisionless interactions are also absent if the 
magnetic fields are weak in the initial fireball. 

Then the relative velocity between pairs and protons 
drives the t wo-stream instability, in particular, of the 
Weibel type (Weibel Il959l) . Let us first consider the lin- 
ear stage of the Weibel instability. The initial distri- 
bution function of each component is given by fj(y) = 

(mj/2TTTj?l 2 e- m i^y- v i )1+ ^l 2T i, where m } is the mass, 7) 
is the temperature, and Vj is the relative velocity for each 
component with j = e + ,e~,p. Initially the system is non- 
relativistic. Linearizing the Vlasov-Maxwell equation, we ob- 
tain the dispersion relation of the pair proton plasma as 

■ ^[1+^-2(0)], (6) 



C 2 k 2 : 



■E 



where ui 2 j = Aimjcfj/ntj is the plasma frequency, Z(£) = 
JJ.-1/2 j°° d x (x— ([)~ l e~ x is the plasma dispersion function, 



and ij = ui/k(2Tj/mj) l/2 dDavidson et all 1 1972b . Note that 
we are considering the Lorentz frame in which the linearized 
Vlasov-Maxwell equation is block-diagonal {D xy = 0), 



^ OJ 2 

j 



Tj k 



(7) 



In the initial stage the relative velocity is small, nijV 2 <C 
Tj, so that we may simplify equation © since ^ < 1. With 
Z(0- ^forf <1, we find to leading order 



1/2 to 2 '-^L-r 2 k 2 



(8) 



Therefore the plasma is unstable (Sw > 0) for < k 2 < k^ 



where k\ ■■ 
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(9) 



(10) 



Equation which characterizes our Lorentz frame is also 
reduced to 



<Y VjOC ^T Vj = - 



(11) 



j 1 j 

Let us apply equation © to the fireball photosphere. Not- 
ing that (a) the pair density dominates the proton density, 
n+ ^n-^>n p (§|2]l, (b) the relative velocity is initially smaller 
than the thermal one, m p V 2 < T+ ~ 7L. ~ T p , and (c) we have 
— n+V + ~ rt-V- ~ ripVp with equation (fTTT i and the charge neu- 
trality J2j qjijVj = 0, we find £\ u^m-Vj/Tj ~ uj 2 pp m p V 2 /T 
and Ej^C^y'C 2 / 7 }) 172 ~ 2uj 2 pe (m e c 2 /T) x l 2 . Then the maxi- 
mum growth rate in equation © is reduced to 

1/2 
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(12) 



where 7Lj = T/O.lkeV and we are assuming V p < (T /m p ) l l 2 . 
Note that the proton current direction is different from the pair 
one, and |V+| ~ |V_| < \V P \ in the frame ( fTTb . 

Comparing the growth rate (3tt>) ma x in equation (fl2l) with 
the acceleration rate ~ g/Vp in equation (O, we find that 
the latter is larger (9w) max < g/V p under the condition V p < 

(T /m p ) A l 2 ~ 1 X 10 7 r_y 2 cm/s. This means that the relative 
velocity V p overtakes the thermal one ~ (T/m p yl 2 before the 
Weibel instability develops. 

Eventually the Weibel instability saturates in the velocity 
region (T/m p ) l l 2 < V p < c. Although no non-linear simula- 
tion exists so far for the proton streaming in pairs, the pro- 
ton and pair flows will split into current filaments (cylindri- 
cal beam) that is bound by the self-generated magnetic fields 
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FIG. 1. — The schematic picture of the unstable GRB photosphere. Ra- 
diation selectively pushes pairs rather than protons, driving the two-stream 
instability between pairs and protons. This produces the small-scale inhomo- 
geneity of the proton-to-pair ratio, that grows up to proton clumps causing 
shocks with pairs. In the accelerating pair frame, particles feel an effective 
gravity with the acceleration g in equation (5). Proton clumps fall down and 
merge, increasing their radii R, infall velocities V p and vertical separation d. 

(ISilva et al.H200l iNishikawa et aUl20Ql lKatoll2005h . Since 
currents with the same (different) direction attract (expel) 
each other, protons effectively separate from pairs 5 and tem- 
porarily produce the small-scale inhomogeneity of the proton- 
to-pair ratio. 

We note that each current can not separate completely 
if the relative velocity V p is smaller than the thermal one 
~ (T /nip) 1 ' 2 because the filament radius ~ fc~' ax in equation 
( [Tol l is larger than the plasma skin depth, r ~ k^ ax > c/w pp , 
and hence the current / = irr 2 qn p V p exceeds the Alfven cur- 
rent I a = m p c 2 V p /q (the maxi mum current limited by the 
self-generated magnetic field) (lArfven|[l939l) . i.e., filaments 
must overlap to reduce th e net curren t. This saturation is 
so-called the Alfven limit (Kato 2005). This is not our case 
since r < c/ u> nn for V p > (T /m p ) x l 2 , i.e., in the particle limit 
dKatol [2005b . Note also that the late evolution is mainly 
driven by protons since the kinetic energy of protons domi- 
nates that of pairs. The generated magnetic fields are weak 
B 2 /8tt < n p m p V 2 < aT 4 at this stage. 

4. POSSIBLE SCENARIO FOR NON-THERMALIZATION 

We have shown that protons will temporarily separate from 
pairs on the small-scale possibly via the Weibel instability. 
The subsequent scenario is very speculative but have attrac- 
tive features that are worth pursuing. We consider in the frame 
comoving with pairs (that is almost the same as the frame in 
equation ^}). In this frame, particles feel an effective grav- 
ity with the acceleration g in equation (0. Unlike pairs that 
interact with radiation, protons fall down, which may lead to 
the "proton sedimentation" in the pair fireball (see Fig . [TJ . 

Initially protons are trapped by magnetic fields on the scale 
of the filament radius. We model the inhomogeneity by the 
proton clumps with the filament size 6 (R ~ cjuj pp ). Since a 
clump consists of many particles, the center-of-mass thermal 
motion is much less than the infall velocity. Also the accel- 
eration time is comparable to the growth time of the instabil- 
ity, so that the clumps continue falling down and collide with 

other clumps within t„, ~ Rj 5V P 



2 x lO^re^V^ s where 
SV p ~ V p is the relative velocity of clumps. Since the Larmor 
radius is comparable to the clump size, the fluid approxima- 
tion begins to be valid. Then, we may apply the Stokes' law to 

5 A small fraction of electrons that preserve the charge neutrality exists in 
the proton filament. 

6 Magnetic clumps are observ ed in the particle simulations, thou gh not 
applicable to our pair proton case. I Chang, Spitkovsky & Arons 2007). 



estimate the terminal infall velocity as V p oc n p R 2 . Since the 
number density of clumps is N c oc n~ l R~ 3 (mass conservation) 
and the cross section for the collision is a c oc R 2 , the merger 
timescale is 

t m ~ (N c <T c 6V p r l oc/T 1 . (13) 
This leads to runaway growth of clumps because the merger 
timescale is shorter for larger clumps. 

After the terminal velocity exceeds the sonic speed of pairs 
Vp > c ( ±, we can not use the Stokes' law but instead the mo- 
mentum balance. Since the pair mass colliding with a proton 
clump during At is AM ~ m e n±R 2 V p At, the momentum con- 
servation gives the velocity change AV P ~ V p AM /m p n p R 3 . 
Equating AV p /At with the acceleration g, we find the termi- 
nal velocity as 

1/2 



v„ 



m p n p gR 
m e n± 



(14) 



With g ~ c 2 /£ p in equation (0, the terminal velocity reaches 
the light speed V p ~ c when the ratio of clumps to the sys- 
tem size is Rjl p ~ m e n±/m p n p < 1. So the final separation 
between clumps in the vertical direction is about 



d~^R 
n 



■I. 



n±m e 



p n p m p 
1a™ ~ , p3. 



(15) 



since the swept mass is R dn p ~ R n p . 

Collisionless shocks arise as the relative velocity between 
proton clumps and pairs approaches the light speed. Pairs 
are shock heated and accelerated to a power-law distribution 
N( r y±)dj± oc n f±d^± for the pair Lorentz factor j± > j m ~ 1, 
where j m is the minimum Lorentz factor. Shocked pairs are 
pushed sideways and then again shocked by other clumps (see 
Fig-Q}- The timescale between shocks is about ~ d/c in equa- 
tion HSi . and interestingly this is comparable to the IC cool- 
ing timescale, 



'y±m e c c n p m p c 



(16) 



d/c cajU-yj^d l±U-y 

for 7± ~ 1 if the proton energy is comparable to the radiation 
energy h p m p c 2 ~ t/ 7 = aT 4 (as our parameters), where we 
use r ie ~ n±OTl P ~ 1 in the second equality. Therefore the 
clumps automatically have the right separation for the contin- 
uous heating of pairs. Since a photon is scattered by j± ~ 1 
pairs once on average for r 7C ~ 1, scattered photons have the 
energy density 

IZUy ~ h p m p c 2 . (17) 
Thus, the non-thermal component has comparable energy to 
the proton kinetic energy, consistent with observations. 

Neglecting multiple scatterings si nce r^ e ~ 1 , we m ay cal- 
culate the observed IC spectrum as (Sari & Esin 200 ll) 

(^/T bs)F v . mdx , V < T b s , 



(18) 



that also peaks at the thermal peak v ~ T y s for p > 2 because 
7,„ ~ 1 (see Figure [2]). The ratio of the IC peak flux to the 
thermal one is ^V,max/-P^max ~ n p w p c 2 /t/ 7 (~ 1 for our pa- 
rameters). At high frequencies v > T t, s , the IC spectrum has 
a cooling spectrum oc v~ p l 2 since the cooling is faster than the 
heating, t c <d/c, for j± > 1. At low frequencies v < T b s , the 
Rayleigh-Jeans spectrum oc v 2 is modified to oc v by down- 
scatterings. Therefore the total (IC plus thermal) spectrum 
is roughly a broken power-law peaking at v ~ T i, s and re- 
sembles the observed GRB spectrum although the low energy 
photon index is slightly harder. 
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FIG. 2. — The schematic picture of the GRB spectrum in our unstable pho- 
tosphere model. Shocked pairs scatter thermal photons into the observed bro- 
ken power-law spectra with energy comparable to the proton kinetic energy, 
consistent with observations. The pair annihilation line is predicted above 
continua in the pair photosphere model, which is blueshifted to v ~ 0.5T3 
GeV and broadened by the order-of-unity distribution of the Lorentz factor 
on the photosphere. At higher energy, a spectral cutoff arises due to the pair 
creation . A closure relation ex ists between observable quantities of lines and 
cutoffs (Murase & Ioka 2007), which may be verified by GLAST. 

5. PREDICTION OF E ± ANNIHILATION LINES 

We can prove that the pair photosphere accompanies the 
pair annihilation line that is detectable above the power-law 
continuum if 7± ~ 1. Let us assume that the continuum ex- 
ceeds the annihilation line, n 7 , e > n 7i /, where n 7 , c (« 7 ,/) is the 
number density of continuum (line) photons at the line energy. 
On the pair photosphere, n±ojl ~ 1, we have n 7 ./ ~ n± be- 
cause the pairs annihilate once o n average in the dynamical 
time, ii±l p /c ~ —n^artp ~ — n± (ISvenssonlll987l) . Then the 
optical depth to the pair production is larger than unity, 

t 77 ~ TUf^arip > n-yjarlp ~ n±ar( P ~ 1, (19) 

and hence the continuum photons create more pairs than the 
existing pairs. However this contradicts with the definition of 
the pair photosphere. Therefore the assumption « 7jC > n 7i ; is 
wrong and the annihilation line should stand out of the con- 
tinuum. As is clear from the above proof, the line flux is com- 
parable to the continuum one when pairs are created from the 
continuum photons. 

The fines are blu e shifte d to v ~ 0.5T3 GeV 
dPe'er. Meszaro s~& Reesl 120061) and broadened by the 
order-of-unity distribution of the Lorentz factor on the photo- 
sphere. At higher energy, a spectral cutoff arises due to the 



pair creation (lAsano & Inoudl2007t iBaring & Hardindll 997 ; 
iLithwick & Sarill200U iRazzaque. Meszaros & Zhangl 120041) . 
A closure relation exists between observab le quantities of 
lines and cutoffs for the pair photosphere ( Murase & Ioka 
12007ft . which may be verified by GLAST. 

6. DISCUSSIONS 

The Yonetoku relation can be derived from equation (fl~|i 
if the photosphere appears at the stellar radius r ~ const., as 
suggested in § [2j and the outflow rate of slow mass satisfies 
M s oc r, because the Lorentz factor after two mass collision is 
about T - (LT S /M s c 2 )^ 2 oc L l l 2 . We can test T oc L 1 / 2 since 
F is determined by the annihilation lines. 

Our model predicts (a) (non-thermal energy) ~ (proton en- 
ergy) ~ (afterglow energy) and hence thermal peaks may out- 
stand for GRBs with > 90% efficiencies, (b) Low energy 
index that is comparable to the observed hardest one. Su- 
perpositions of pulses or incomplete thermalization may pro- 
duce softer ind ex, (c) Dim optical flashes from reverse shocks 
dLi et al.ll2.Q03b . consistent with recent observations, (d) Po- 
larization in the non-thermal component (except for the ther- 
mal peak) produced by IC scatterings, (e) High energy cos- 
mic rays and neutrinos since protons might be also acceler- 
ated at shocks. Particle simulations with pairs and protons are 
needed to verify our model. Other plasma instabilities, such 
as the electrostatic two-stream instability or the drift-wave in- 
stability with the initial magnetic fields, may also generate the 
initial proton inhomogeneity. 



We thank K. Asano and T. Kato for useful comments. This 
work is supported in part by Grant-in- Aid for the 21st Cen- 
tury COE "Center for Diversity and Universality in Physics" 
from the Ministry of Education, Culture, Sports, Science and 
Technology (MEXT) of Japan and by the Grant-in-Aid from 
the Ministry of Education, Culture, Sports, Science and Tech- 
nology (MEXT) of Japan, No. 18740147 (K.I.), No.19104006, 
No. 19740139 (S.N.) and No. 19540283, No. 19047004 (T.N.). 



REFERENCES 



Alfven, H. 1939. Phys. Rev., 55, 425 

Amati, L. 2006, MNRAS, 372, 233 

Asano, K., and Inoue, S. 2007, arXiv:0705.2910 

Baring, M. G., & Harding, A. K. 1997, ApJ, 491, 663 

Beloborodov, A. M. 2000, ApJ, 539, L25 

Chang, P., Spitkovsky, A., and Arons, J. 2007, arXiv:0704.3832 

Cox, A. N. 2000, Allen's Astrophysical Quantities (4th ed.; New York: AIP) 

Davidson, R. C, Hammer, D. A., Haber, I„ and Wagner, C. E. 1972, 

Phys. Fluids, 15, 317 
Ghisellini, G., and Celotti, A. 1999, ApJ, 511, L93 
Ghisellini, G., Celotti, A., and Lazzati, D. 2000. MNRAS, 313, LI 
Giannios, D, and Spruit,H. C. 2007, A&A, 469, 1 
Goodman, J. 1986, ApJ, 308, L47 

Ioka, K., Toma, K., Yamazaki, R., and Nakamura, T. 2006, A&A, 458, 7 

Kato, T. N. 2005, Phys. Plasmas., 12, 080705 

Kobayashi, S., and Sari, R. 2001, ApJ, 551, 934 

Li, Z., Dai, Z. G., Lu, T, and Song, L. M. 2003, ApJ, 599, 380 

Lithwick, Y, & Sari, R. 2001, ApJ, 555, 540 

Meszaros, P. 2006, Rep. Prog. Phys., 69, 2259 

Meszaros, P., Ramirez-Ruiz, E„ Rees, M. J., and Zhang, B. 2002, ApJ, 578, 
812 

Meszaros, P., and Rees, M. J. 2000, ApJ, 530, 292 
Murase, K., and Ioka, K. 2007, submitted 



Nishikawa, K.-I., et al. 2003, ApJ, 595, 555 
Paczyriski, B. 1986, ApJ, 308, L43 

Pe'er, A., Meszaros, P., & Rees, M. J. 2006, ApJ, 642, 995 

Pe'er, A., and Waxman, E. 2004, ApJ, 613, 448 

Razzaque, S., Meszaros, P., and Zhang., B. 2004, ApJ, 613, 1072 

Rees, M. J., and Meszaros, P. 2005, ApJ, 628, 847 

Ryde, F. 2005, ApJ, 625, L95 

Sari, R., and Esin, A. A. 2001, ApJ, 548, 787 

Silva, L. O., et al. 2003, ApJ, 596, L121 

Svensson, R. 1987, MNRAS, 227, 403 

Thompson, C. 1994, MNRAS, 270, 480 

Thompson, C, Meszaros, P., and Rees, M. J. 2007, ApJ, in press, 
astro-ph/0608282 

Toma, K., Ioka, K., Yamazaki, R., and Nakamura, T. 2006, ApJ, 640, L139 
Waxman, E., and Piran, T. 1994, ApJ, 433, L85 
Weibel, E. S. 1959, Phys. Rev. Lett., 2, 83 

Yonetoku, D., Murakami, T., Nakamura, T, Yamazaki, R., Inoue, A., and 

Ioka, K. 2004, ApJ, 609, 935 
Zhang, B. 2007, Chinese J. Astron. Astrophys., 7, 1 
Zhang, B., et al. 2007, ApJ, 655, 989 

Zhang, W., Woosley, S. E., and Heger, A. 2004, ApJ, 608, 365 



